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We consider some aspects of the global evolution problem of Hamiltonian homogeneous, anisotropic cosmologies 
derived from a purely quadratic action functional of the scalar curvature. We show that models can isotropize in the 
positive asymptotic direction and that quadratic diagonal Bianchi IX models do not recollapse and may be regular 
initially. Although the global existence and isotropization results we prove hold quite generally, they are applied 
to specific Bianchi models in an attempt to describe how certain dynamical properties uncommon to the general 
relativity case, become generic features of these quadratic universes. The question of integrability of the models 
is also considered. Our results point to the fact that the more general models are not integrable in the sense of 
Painleve and for the Bianchi IX case this may be connected to the validity of a BKL oscillatory picture on approach 
to the singularity in sharp contrast with other higher order gravity theories that contain an Einstein term and show 
a monotonic evolution towards the initial singularity. 



1. Introduction 

The problem of describing the past and future asymp- 
totic states of general cosmological models occupies 
a central position in mathematical cosmology. This 
question can be addressed either in the framework of 
general relativity and string theory or in that of higher- 
order and scalar-tensor theories. Among several dif- 
ferent approaches that have been developed over the 
years to tackle the asymptotic problem, two are dis- 
tinguished, namely, the Hamiltonian approach and the 
dynamical systems approach. Both have been and are 
being used widely and with great success when dealing 
with the asymptotic problem in the framework of gen- 
eral relativistic and string cosmology but less widely 
so in the case of higher-order cosmologies, that is, in 
the description of the asymptotic states of cosmolog- 
ical models derived from an action functional which 
depends nonlinearly on the curvature invariants. 

In this paper we use the Hamiltonian formulation of 
homogeneous but anisotropic cosmologies derived from 
an action integral purely quadratic in the scalar cur- 
vature (quadratic cosmologies) developed by Demaret 
and Querella in [|j to prove several results pertaining 
to the global structure of the solution space of these 
systems. The main motivation for studying such sys- 
tems is not so much due to their confrontation with 
the currently available empirical data (see, however, 
||) but to their different nature. 

It is certainly true that recent high-precision tests 
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such as the Hulse- Taylor binary pulsar (cf. have 
confirmed general relativity to an outstanding degree. 
This theory predicts that in the case of the binary 
pulsar PSR1913+16 the orbit slowly decays and the 
orbital period increases due to the loss of energy by 
emission of gravity waves and this has been observed 
with high precision. However, general relativity leads 
to singularities in the spacetimes of all known cos- 
mological models - places where further predictions 
about the structure and evolution of the universe can- 
not be made. Quadratic curvature terms present in the 
gravitational action, which naturally arise in string- 
theoretic considerations as quantum gravitational cor- 
rections, may rectify this situation and lead to cos- 
mological models free from such pathologies. One 
then envisions the general-relativistic (linear) action 
to be a low-energy, low spacetime curvature limit of 
some quantum gravity theory that would contain such 
higher-order curvature invariants and which has yet to 
be worked out. Of course the purely quadratic theories 
considered here do not reproduce the weak field limit 
of general relativity. However, such theories may gain 
importance in the high-energy regimes of the early 
universe and serve as prototypes from which more 
elaborate or natural, Lagrangian theories of gravity 
can be built and studied. 

Due to their conformal equivalence [JjJ , the dynam- 
ics of higher-order vacuum cosmologies can be regarded 
as closely related to that of general scalar field cos- 
mologies in general relativity. However, it should be 
emphasized that this procedure, although practically 
useful, cannot imply a substitute for a full analysis 
of the original (higher-order) system since, for exam- 
ple, at those points on the manifold where the con- 
formal transformation is singular, one cannot obtain a 
complete picture of the original dynamics when invert- 
ing the conformally related one. Therefore one has to 
rely solely on the original (higher order) field equations 
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when analysing a higher-order system. 

We consider vacuum Bianchi cosmologies of class A 
in the expanding direction and prove global existence 
results towards a future asymptotic state. We show 
how these quadratic cosmological systems isotropize 
and prove that Bianchi IX models in the expanding 
direction have a monotonically increasing volume and 
hence cannot recollapse. This result is in sharp con- 
trast with the known behaviour of these systems in 
general relativity (cf. 10 , UW). The question of the 



algebraic integrability of quadratic cosmologies is also 
analysed by applying the so-called Painleve test to the 
relevant cosmological equations, and their Lie point 
symmetries are calculated. 

The plan of this paper is as follows. In the next 
section we cast the purely quadratic homogeneous cos- 
mologies in a Hamiltonian form (cf. || ) and write down 
the basic Hamiltonian cosmological equations. This 
Hamiltonian reduction process has an added advan- 
tage to that with the usual Lagrangian field equa- 
tions, namely, the original fourth-order (Lagrangian) 
field equations are reduced to a first order system in 
the Hamiltonian variables. In Sec. 3, by defining new 
variables, we further reduce the Hamiltonian equations 
to a suitable autonomous dynamical system for the vol- 
ume and anisotropy functions. Sec. 4 is the heart of 
the paper. We prove two theorems about the global 
existence of solutions to the basic equations in the ex- 
panding direction. These theorems are then used in 
Sec. 5 to prove an isotropization theorem for a class of 
quadratic, homogeneous, anisotropic cosmologies, and 
all previous results are applied in Sec. 6 to several spe- 
cific Bianchi metrics in an attempt to highlight their 
dynamics. In particular, we show, as a corollary of 
the basic theorems, that there exist non-recollapsing, 
ever-expanding Bianchi IX cosmologies in marked con- 
trast to what is known about such systems in general 
relativity. We round up our analysis in Sec. 7 by list- 
ing the Lie point symmetries of quadratic cosmologies 
and performing for them the Painleve test to decide 
on their algebraic integrability. We conclude in Sec. 8 
wherein we discuss prospects for future work. 

2. Purely quadratic Hamiltonian 
cosmologies 

In this section we derive the general expressions for 
the super-Hamiltonian, supermomenta and canonical 
equations which are our starting point for the analysis 
of relevant cosmological models in subsequent sections. 
Our beginning is the purely quadratic action, 



(1) 



Following the ADM Hamiltonian prescription, we as- 
sume a foliation of spacetime into spacelike hypersur- 
faces and use the ADM coordinate basis to compute 



the components of various curvature tensors. The "0" 
index below refers to the normal component to the hy- 
persurfaces, while the superscripts " and indicate 
quantities defined over the 4-dimensional spacetime 
and the spacelike hypersurfaces, respectively. Most of 
the time, when there is no ambiguity, the last subscript 
will be omitted. This splitting of spacetime allows one 
to rewrite the action ([!]) in the form 



(2) 



with 



(4) P = (3) P + K 2 - 3 tr K 2 - 2 



V 2 N 
N 



-2g kl C n K k 



ki- 



CrKm = -{K klfl - N*K k ni - N) k K u - N^K ik ), 
11 

where K%j and ^'Rij are the extrinsic and intrinsic 
curvature tensors relative to the hypersurfaces men- 
tioned above, the symbol | denotes covariant differen- 
tiation on these 3-surfaces, N and 2Vj are the lapse 
and shift functions and denotes the Lie derivative 
along the normal ft to these hypersurfaces. 

The Hamiltonian form of (Q) is then given, up to 
surface terms, by the integral 



S = 



with 



J d A x N[Pij £„•(/•' • //•' C n9ij -n(g, Q,p, P)](3) 



H(g,Q,p,P) = -2p»P k 



Q 2 



qVr.. - Q R+ QfiPyP + |(trp 2 - p 2 ), 



(4) 



P — K 



(5) 



where the trace and the traceless part of a tensor A$j 
are denoted respectively by A and Afj , tr A 2 — A^Aij 
and 



H(g, Q, K) = KQPKij + R tJ + Q v {ij 



Q 2 , Q 



IzPcg 1 ' 2 2 



(trK 2 - K 2 - R), (6) 



wherein Q lJ is a tensor introduced in order to re- 
duce the order of the Lagrangian as in Ostrogradski's 
method. The canonical variables are gij and Q lJ and 
their conjugate momenta are p lJ and Py , respectively. 
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We can also write the action (||) in the usual Hamil- 
tonian form, where all the constraints are manifest, 
namely, 



S = 



P^Q 13 + P l] g I3 -n*(g,Q lPl P) , (7) 



where we have introduced the total Hamiltonian con- 
straint H* as 



H* =NH + N k H k = N^H^ 



(8) 



TL and TL k being the super-Hamiltonian and supermo- 
menta (or spatial constraints) given by (||) and by 



H k = Q ij P ij]k + 2{P ik Q% j - 2g lkP 11 



(9) 



respectively. 

The 3-metric of any diagonal Bianchi-type model is 

da 2 = e 2 »[eW++^-\u 1 ) 2 + eW+-^-\L0 2 ) 2 

+ e" 4 ' 3 +(^ 3 ) 2 ], (10) 

where fi, f3+, /3_ are functions of time only and the set 

{Q 1 } is the basis of 1-forms with structure coefficients 

ni 



dQ % = \C) k uj' j A u 



Following 0, we perform a canonical transformation 
from the original set of canonical variables {<?, Q, p, P} 
to the set {/x, (3 + , /?_; U^, II+, n_; Q + , Q_, Q n , P + , 
P- , P n }, where 



9n 

522 
#33 

n\ 
n 2 2 
n 3 



e 2 Ate 2(/3 + + v / 3/3_)^ 

e 2 AIe 2(/3 + -V3/3-) ; 
e 2M e -4/3 + 



^-(2n Al + n+ + V3n_), 
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-L(2n M + n + -^n_), 



(Hu-n + ), 



P,:, = 



II'' + P T £Q Tjfc 



pT. + —o 



x/3 

, p>q 



x/3 



Q Tl ' 



V3 



p 1 2 _ 



-P diag (Q+ + V3Q_ , Q+ - V3Q_ , -2Q, 
V6 

diag(P + + v^P-, P+ - V3P_, -2P+). 

V6 



In terms of these new variables the original Bianchi 
action can be written as 



S 



d 4 n (n> + n+/?+ + n_/?_ + p+q+ 



n* 



+P_Q_+P n Q, 

with the 4- volume element d 4 fl = dt Auj 1 Auj 2 Aui 3 . The 
explicit form of the total Hamiltonian constraint Ti.* 
depends on the particular Bianchi model considered. 
In the case L = R 2 we have to impose the constraint 
Q ± ?a . The super-Hamiltonian then reduces to 



H 



n 2 



6A 



-Qr 



V3 



2 



P n Qn 



(11) 



where V* = V*(0+, is directly related to the 

usual Bianchi potentials (cf. JT^ j . ch. 10). Denoting in 
everything that follows by X y the derivative of the 
function X with respect to the variable y , the canon- 
ical equations obtained from dill) are 



A* = 



#t = 



-^P 

at n± 

2V3 Qn^ 

1 



n " - " 2ft 

AT 

n+ = =e 



-3^n 2 



Qn 



-2// 



QnV* 



V3 



/3±' 



Pn = 



N 

4\/3 



n+ 



"3/3 c e " 3AlQ " 



n 



8Pi 



(12) 



3. Reduction to an autonomous form 

By imposing the weak constraint TLr sa and manip- 
ulating the basic Hamiltonian system ( |l2| ) we obtain 



PnQn ^ H/i 



(13) 



where k is a constant of integration. When the scalar 
curvature is not constant, we can fix the temporal 
gauge by choosing = —t to find 



Qn — 2 fic& 



N 



3& 
2k 



(14) 
(15) 



By some further straightforward algebraic manipu- 
lations it is possible to reduce the set of canonical 
equations to a non-autonomous dynamical system 
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of three coupled equations for the physical variables 
A*(t), 0±(t): 



p!'t + fj,' + A^c 6 ^ - A 2 tc 4 ^V = 0, 
P" t + (3' + A 3 te^Vf} = 0, 

where Ai are strictly positive constants given by 

9p , 3p 3p Pi 



At 



16' 



A 2 



A, 



8 ' 



P 



(16) 



;i7) 



with 1 and fc 2 arbitrary constants. (The canonical 
variables fj, and (3± become uncoupled when the poten- 
tial term vanishes, i.e. for Bianchi-type I model.) We 
can rewrite the dynamical system (|16|) in the so-called 
Taub gauge, d\ = e^^dt, where A is the physical 
time with < A = A(t) < oo and, assuming invert- 
ibility, t = t(X), we can write fi(t) = yu(A) , (3+(t) = 
/3+(A) or 0_(t) = /3_(A) and V = V{fi+(t),0-{t)) = 
V(/3 + (\),P-(X)). 

The treatment of the system of ODEs ( |l6|) is greatly 
facilitated by transforming it into an autonomous sys- 
tem. For this purpose we introduce new variables 
(t* ,w,y,z) as follows: 



0-(t) = z(f), 



0+V)=y(t*), 



(18) 



where tp is a smooth function of t* . On inserting dl8| ) 
into @ we obtain: 



wipip 



~ + A ie 6w <p 2 

-A 2 e 4w V V(y,z) = 0, 

[3<p fitpCp 4 

-rr rr + ~ + ^3e V>(y,«)^ = 0, 



(19) 
(20) 



where the dot denotes differentiation with respect to 
t* and 



= !/(**) 



*(**). 



(21) 



In (|19|) and ( |20| ) we eliminate w and /?, respectively, 
by requiring that 



1 

and hence 



<PjP 



ip{t*) = Ce fcl 



(22) 



(23) 



Eq. ( J23| ) is the general solution of ( |22| ) with C and fci 
arbitrary constants. On the substitution of fl23| ) into 
@ and @ we deduce that 



C/fc 2 



+ AiC 2 e 6w+3klt 
+ A 3 V fj (y,z)Ce 4w+2k ^ 



A 2 V(y,z)Ce 4w+2klt " = 0, 
(24) 



From (E4J) it is obvious that by introducing a new de- 
pendent variable x(t*) through 

w{t*) = ±[x(t*) - kit*] (25) 

and, since it suffices without loss of generality and for 
our purposes, by taking in (p3|)-(p5|) 



C = h = -1, 

the basic system 
form 



(26) 

finally assumes the autonomous 



x + + a 2 e 2x V(y, z) = 0, 

y + a 3 e 2x V y (y,z) = 
z + a 3 e 2x V z (y,z) = 0, 
H(t) = ±[x(f) + f], /3+(t) 
0_(t)=z(t*), t = -< 



y(t* 



(27) 
(28) 
(29) 

(30) 



with a\ — —2Ai , a 2 — ~2A 2 , a 3 = A 3 . All functions 
in this system can also be regarded as functions of A . 

Notation: In the following we use the two basic time 
variables t* and A . A prime will denote differentiation 
with respect to A and an overdot differentiation with 
respect to t* . 

4. Global existence 

In this section we show that under very plausible hy- 
potheses on the potential V(y, z], the solutions of the 
autonomous dynamical system d27j)-(|29|) can be defined 
and smoothly extended in the expanding direction for 
the whole positive A-halfline. The main result of this 
section is given in the following two theorems. The- 
orem 1 basically says that t* is a compactified time 
parameter meaning that its interval of definition is a 
compact interval on the real line and Theorem 2 proves 
that the interval of definition of the time variable A is 
the interval [0, oo) . 

Theorem 1. If the potential V = V(/?+(A),/3_(A)) 
is a smooth function of (3+(X) and /3-(A) and satisfies 
for all real (3+(X) and /3_(A) the condition 



V = V(f3 + (\), /3-(A)) > M, Mel, 



(31) 



and in the case M < the arbitrary constant p > 
given in (plf) fulGls 



P < 



36(4) 



27e 3x ° - 256M 3 ' 



(32) 



then there exist solutions /j(X) and (3 + (X) , /3_(A) to 
([l6|), which are monotonically increasing and decreas- 
ing, respectively, defined on the interval < A < A max , 
where 



A max — 



o.xp <; -x[- ln(-w)] \ (~w)- 3/2 dw, (33) 
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T* is a finite number and xq , x and t$ have been 
introduced in (fl£ 



Theorem 2. If Theorem 1 is valid and, in addi- 
tion, in the case lim,\— >A max /3+(A) = — oo and (or) 
limA^A max /3-(A) = — oo, the derivatives Vp are abso- 



lutely bounded, that is, for n > 1 , fci G M 
G R", /3i m) G R", and /3+(A) G 
/3-(A) G [— oo,/?_ ], tie conditions 



hoid, then necessarily 

Amax OO, 



lim fJ-(X) = oo. 

A — »-oo 



• , fc 2 G R+ 
-oo,^ m) ] 



(34) 
(35) 

(36) 



We break the proof of these theorems into the two 
Lemmata below. We define 



Pi =X, p 2 = V, P3 = z. 

Then @-(||) become 

PiPi,x + aie 3x + a 2 e 2x V{y, z) = 0, 
P2P2, y + a 3 e 2x V y = 0, 
P3P3,z + a 3 e 2x V z = 0, 



(37) 



(38) 



and after integration the system ( |3q ) yields by virtue 
x = (^A ie 3x +4A 2 J e 2x V(y,z)dx^j ' , (39) 

r \ 1/2 

-2A 3 J e 2x V y dy\ , (40) 

1/2 

(41) 



= - [ -2A 3 I e 2x V r dz 



Each choice of positive and negative signs in front of 
the square roots in (|39|) and (f!o[) is taken in order 
to generate increasing x(t*) and decreasing y(t*) and 
z(t*) as functions of t* , respectively, and will be justi- 
fied at the end of this section. Since the system (27)- 
( p9| ) is autonomous, the initial value of t* , ij*, can be 
arbitrary, and we consider here t* = tn. Therefore, 
by incorporating into our system (p7|)— (29) the physi- 



cally plausible assumption that V(y, z), V y (y,z) and 
V z (y,z) are continuous for all y and z in the vicinity 
of t* = 4q , we can invoke the existence theorem for 
the dynamical system (|27|)-(p9|) and establish the va- 
lidity in a neighbourhood TV* of t* = t$ of a unique, 
C 2 solution (x(t*),y(t*),z(t*)), 

x = x(t* ,x ,y ,z ,x ,2/o,Zo), 
y = y(t*,xo,yo,z a ,xo,yo,z ), 

z = z(t*, x , y , z 0) x , y ,z a ), (42) 



satisfying the initial conditions 

x(tl) = x , x(tj$) = x > 0, 

y(*o) = vo, y(to) = m < o, 

z (*o) = z o, «(*o) = i o < 0, 



(43) 



where the numbers Xq, , yo, z , are arbitrary and Xq, yo , 
z are signed as shown, x(t*) and y(t*) , z(t*) being 
monotonically increasing and decreasing functions of 
t* , respectively. 

Equivalently, due to the monotonicity of x(t*), we 
may choose x as the independent variable and express 
the solution (W2h as 



y = y{x), 



z(x), 



(44) 



with the roles of t* and x interchanged in the first of 
Eqs. (fi^). Analogous to relations can be written 
down by considering y or z as independent variables. 

Before we consider the existence proof, we proceed 
to compactify the time interval [0, A max ] . From (|3£~ 
( fill ) we deduce by using © and (H) 



T* - t* = 



Xq < x(t*) < X\ = x(T*) 



(45) 



where 



6p I e 2w V(y(w),z(w))dw, 
yi 



T*-t* = - 



H 



-1/2 



yi = y(T*) < y(t*) < y ; 

yo 

here H(y) = y 2 + (3/4)p J e 2x ^V y dw and 

y 



(46) 



T*-t*=- 



zi = z{T*) < z(t*) < z , 



(47) 



where K(z) = z 2 + (3/A)pje 2x< - w ^V z dw. In @ the 

z 

upper limit x\ in the integral denotes the maximum 
x = x(t*) — the value for which the solution ( p4| ) is 
defined, i.e., xq < x < x\ , xq < w < X\ , which implies 
that the solution (j5|) holds for t% < t* < T* with x x = 
x(T*), yi = y(T*j = y(xi) and z\ = z(T*) = z(xi) . 

We observe now that by virtue of dA = e 3p ( t Mt 
and © 



A(t) = 



t=-e _ 



(48) 



and that A max is given precisely by the form in (p3[). 
Therefore, by means of (|30| ) and (|48|) we have defined 
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a bijective mapping of the physical interval < A < 
Amax onto fg < t* < T* . For this mapping to be 
physically legitimate and useful for the treatment of 
( |l6| ) we must first prove (for the class of potentials 
V(y, z) of interest entering (jl^) or, equivalently, ( |2^ ) 
(|g)) that 



Amax — 



(49) 



In conjunction with ( j49| ) we observe that the inte- 
grand in (H|) can diverge only at — e _T * . This occurs 
either for lim t *^T* x(t*) — oo, irrespective of whether 
T* is finite of infinite, or for T* = oo (w = 0). In both 
cases (Q holds. 

We first consider the two possibilities regarding the 
values of T* , i.e., finite or infinite, and in this respect 
state the following lemma for (|45|)-(^7|), which is in 
fact valid for a class of physically interesting potentials 
V(y,z). 

Lemma 1. If the potential V = V(/3 + (X),/3-(X}) 
satisfies for all real (3+(X) and /3_(A) the condition 



V = V(f3+(\),(3-(\))>M, 



m e 



(50) 



and in the case M < the arbitrary constant p > 
given in (flTf) fulfils 



P < 



36 (x' o y 



27e 3x ° - 256M 3 : 



(51) 



xq and x having been introduced in (j43|) ; then T* 
and thereby the maximal interval of existence, [£q , 
T*] , of solution are finite. 



Proof. Since by assumption V > M , Me R, for 
all real /3 + (A),/3_(A), then due to (§0|) also V(y,z) > 
M , M G R, for all real y and z, and we obtain the 
estimate: 



G(x) > (xo) 2 



(io) 2 



1 PB 



•.i h 



3pMe 2xo 



6pM 



e 2w dw 



-pe 2x (e x +4M), 



(52) 



and we call the expression appearing in the last line 
F(x) . Now to further simplify the notation, let us call 
F Xo (x) the 'value' of F(x) when the exponential in the 
fourth term multiplying the round bracket is e 2x ° . We 
distinguish the following cases: 
A) M + e* /4 > 0. 

Then, since x = x(t*) is an increasing function of 
t* , the inequality M + e x /4 > is valid and, owing to 
(j52|), for x > x 



F > Ft 



( x q) 



-pe 3x ° 
4 



-pe 2xo e x > 0. 
4 



(53) 



Eqs. (Kg) and <M) yield for (|5l) the estimate 



< T* - t* 



G -i/2< / F -y 



< / F-^ 2 =J Xo {x x ). 



(54) 



B) M + e xo /4 < (M < 0). 

Consider /(x) = | pe 2:c (4M + e^) . 
Bi) -8M/3 < e X0 < -AM. Then /(x) increases for 
x > xq . Thus 

minF = F(x ) = (x' ) 2 > 0. 

Therefore 

F(x) > 0, X > Xq 

and due to ( |52]) we obtain for ( ^5| ) 

X! ln(-4M) 



< T* - t* = G< 



where 



! t j / Fd.T + Ji(xi), (55) 

X 2!0 



Ji(a;i) 



^ln(-4M)( a; l)- 



(56) 



B 2 ) e 210 < -8M/3. Then f(x) = \pe 2x (AM + e x ) 
decreases for xq < x < ln(— 8M/3) and attains its 
minimum value for x = ln(— 8M/3), i.e., min/(2:) = 
64pAf 3 /9 < 0. Consequently, 

minF = (± ) 2 - 3pMe 2xo - lpe 3x ° + fpM 3 > (57) 

by choosing the arbitrary constant p > in ( p7| ) so 
that ( pl| ) holds. Thus we ensure that F(x) > 0, x > 
xq and yet again arrive at (|56|). Since the integrals 
J Xo (oo) and Ji(oo) in ( |54| ) and (p6|), respectively, exist 
H p. 92, 2.315], the lemma is proved. □ 

We note that by construction (f|5[)-(|47]) as well as 
the estimates (|^) and (|55| ) hold for all xo , xo > , yo , 
yo < and zo, io < 0. Also the condition ( fill ) is 
not particularly restrictive since it involves a relation 
between the arbitrary constants p > 0, and io > 0. 

Now, by virtue of Lemma 1 and the discussion fol- 
lowing ([49|), we observe that the remaining possibility 
for the validity of ( f49| ) is that lim t »^x* x(t*) = 00. 
Consequently the ensuing developments focus on the 
proof of this fact for all xq, xq > , yo, yo < and 
z , i < . 

To this end we use an immediate and important 
consequence of Lemma 1, namely that |j) at least one 
of lim t ,^r. |a;(f*)|, lim t ,^r» \y(t*)\, lim t .^ T * \z{t*)\ 
is 00 . If for a given V{y, z) the analytic structure 
of either one or of both of the quantities dV(y, z)/dy 
and dV(y, z)/dz is such as to exclude in the context of 
( p6[ ) and (|jj) for all initial conditions x , io > , yo, 
yo < and z , z < the cases y = —00 or z = —00 
or ?/ = —00 and z = —00, then clearly in ( ^6|) and 
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(eZD 2/i ^ R an< i zi whence necessarily lim t »^T* 
x(t*) = oo . 

If, however, the previous cases cannot be as above 
a priori eliminated, upon taking into account the an- 
alytic form of some physically interesting potentials 
V(y, z) , we have the following 

Lemma 2. Suppose that for the initial conditions 
xq, xq > , yo, yo < and zq, zq < , so that either 
yi = — oo or z\ — —oo or y = — oo and z = — oo 
appear for either w G [— oo,y m ] or w € [—oo,z m ], 
where y m G R~ and z m G M~ , the derivatives V y , V z 
in (p^j and ( fi^ ) are absolutely hounded, that is, for 
G M+,n > 1 



|V>"| < fci 
and 

\V z w n \ < k 2 . 



(58) 



(59) 



Then, for all initial conditions xq , xq > , yo , yo < 
and zq, zq < we have 



lim x(t*) = oo. 

t *_>T* 



(60) 



Proof. The proof is carried out by reductio ad absur- 
dum. Assume firstly the existence of initial conditions 
xq, xq > , yo , yo < and zq , Zq < so that con- 
currently 

lim x(w) = lim x(t*) = x\ G R, 



lim z(t*) = zi G 



and 



lim y(i* 



2/1 



-oo. 



We use @ and @ to obtain 
Hm^ = -ff 1/2 (-oo). 



(61) 
(62) 



Obviously now for y?n < 0, owing to ( |5q ) and since 
xq < x(w) < x\ for — oo < w < yo j we obtain for the 
integral in (pa) 



2a;(u;)y 



< 



!/() 



3 2an(to)^ 



(/((> 







rvo 


e 2x{w) V y 


oo a "* 


dw + 











dw 



< 



fix{w)y 



dw + kie 2xi 



l dw. 



(63) 



Due to the continuity of V y as a function of y = y(t*) 
and z — z(t*) = z{y) (owing to the monotonicity of 
y(t*)), with -oo < y(t*) < y and z x < z(t*) < y Q , 
the integral 



rvo 


e 2x{w) V y 


dw 


>y m 







exists. Consequently, the integral on the left-hand 
side of (|63|), by virtue of the estimate on the right- 
hand side of (^|), also converges. Therefore ( |62|) yields 
limt*_>T*2/ = —ax, 0,1 € R + , ai depending of course 
on xq, xq > 0, yo, yo < 0, zq, zq < 0. This relation im- 
plies that in some left neighbourhood of T* we have 

y = -a(t*,x ,x > 0, yo,yo < 0, z ,z < 0) < 0, 

with \a(t* ,xq,±q > 0,j/0) yo < 0,zq,z < 0)| < oo and 
limt*—,^. a(t*) = a\ , whereupon by integration we de- 
duce that in fact lim t *^r* y(t*) is hnite. 

We thus arrive at a contradiction with our hypoth- 
esis for the existence of appropriate xq , xq > , yo , 
2/o < and zq, zq < such that Unx u) ^ J/1= _ 0O x(w) — 
lirn t *_> r * x(t*) = xi G R, lim t »^. r » = Z\ G R 
and ]im t ._ >r . = 2/1 = -00. 

Likewise, due to (|59|), it is verified that there do 
not exist Xo , xq > , 2/0 > Vo < an d 2 , io < such 
that liimjj^^^-oo a;(u>) = Mva.f'^T* x{t*) = x\ G R, 
lim t *^T* = 2/i and lim t »^T* = zi = —00 

or, owing to ( p8| ) and (p9[), such that limm^-oo x(u>) = 
lim t »^ r » x(t*) = .xi G M, lim t »^ T » = 2/1 = -00 
and lim t *^r» z(t*) = —00. 

These results, according to the remarks prior to 
Lemma 2, immediately imply (|60|). The proof of 
Lemma 2 is now complete and with this, the proofs of 
Theorems 1 and 2 are also complete. □ 

At this stage some remarks are appropriate. The 
conditions (^) and ( [59"1) are sufficient only, as we shall 
later demonstrate, in the case of the Bianchi IX po- 
tential. Secondly, the validity of ((6^) does not either 
entail limt»^T* y(t*) = —00 or limi»_>T* z(t*) = —00 
or both, or preclude these possibilities, which may 
appear for certain values of Xo , xq > , yo, yo < 
and zq , zq < or for specific analytic forms of 
V(y,z) (cf. also the proof of Theorem 3). Essen- 
tially it has been verified that lim{»_>r* y(t*) = —00 
or lim t »^x* z(t*) = —00 or both are not possible with- 
out the validity of lim{»_>r» x(t*) — 00. 

On recalling the discussion following (^9|) we de- 
duce that Lemmas 1 and 2 verify ( |49| ) and, therefore, at 
least for the class of potentials V = V((3+(X), (3 (A)) = 
V(y,z) satisfying conditions (p0|), fl5S| ) and ( |5^ ) we 
have in fact mapped < A < A max = 00 onto t^ < 
t* < T* , T* being finite. In this sense we have com- 
pactified the interval [0, 00) in which A is defined. 

To sum up, by utilising the transformations (|3 
and dA = c z ^At, we have 



(j,' = i(i + l)e t *- 3 ^ l(A) >0, 
+ = 2/e t *exp{-|[x(r)+r]} < 0, 
0'_ = ie**exp{-|[a;(r)+r]} < 0. 



(64) 
(65) 
(66) 



Consequently, by virtue of (|3C|), (|3J)-(|41|) and (|3|, 
©-© and since y{t*) and z(t*) exist for % < t* < 
T* , we conclude that the functions ;u(A) and /3±(A) 
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are defined for < A < oo and that they are mono- 
tonically increasing and decreasing functions of A , re- 
spectively, fulfilling the initial conditions 

50 c °+*o)> 

l(x + 1) exp[t* - |(x + t*)} > 0, 



/to 
Mo 

P+,o 



2/o, 

2/o exp[t„ 



i exp[^ - |(x + to)] < 0, 
whereby, also by virtue of @), (||), ([l' 
oo, Mo < m(A) < oo, 



(67) 
(68) 



(69) 



and M), 



lim fi(X) 

A — >oo 



r* - 1* = 



lim x(i* 



xi = oo; 



(70) 



lim /3+(A) = j/i, 2/i < /3+(A) < y , 



2/i ; 



(71) 



lim /3_(A) = #i 5 

A — *oo 



zi < /3_(A) < z , 



T* - 1* = 



lim z(T*) 



K 



-1/2 



(72) 



We stress that (|67() - (|72|) remain valid for aiZ xo , ±o > 

, yo , 2/o < and z , z < . 

5. Isotropization theorem for 
quadratic cosmologies 

In view of the results of the previous section, it is now 
possible to prove an isotropization theorem for those 
cosmologies that have potentials complying with the 
assumptions of Theorems 1 and 2. We consider the 
scalar shear parameter £ and show that it becomes 
vanishingly small in the expanding direction. We have 

Theorem 3. The validity of Theorems 1 and 2 im- 
plies that the shear fulfils 

lim £ 2 (A) = lim U/3' + /^') 2 + (f3'_/n') 2 ] = 0. (73) 

A — >oc A — *oo 



Proof. We first observe that due to (64), fl65| ) and 

(74) 
(75) 



/?;/// = 2y/(x+l), 
&_lvt = 2z/(x + l), 

with (cf. di), @) 

x = G 1 / 2 . 



(76) 



On recalling that liim^^ e 2x ^ ~ oo , we consider 
( [iii"l) and observe that the following possibilities exist: 



1) The integral in (|6l]) for w — > y\ diverges, in which 
case (cf. ( |58| ) for the case 1/1 = -oo or zi = -oo 
or yi = — oo and zi = — oo) at the most 



lim y 



lim e x{t ">. 



(77) 



2) The integral in (|61|) for w — » y\ converges and thus 



lim y = —a. 



a € 



(78) 



Let it be pointed out that by using standard crite- 
ria for the convergence of integrals one can easily con- 
struct examples for possibilities 1) and 2) above. Cases 
( \n\ ) and (|7q ) , which are of course compatible with (|5^ ) 
(i.e., (|34|)), are the ones possible for the behaviour of 
lim t »^T* y, in view of y(t*) being monotonically de- 
creasing for t* G [ £q,T*), with (77) corresponding to 
lim t *^T* y(t*) = — oo or lim t »^T* y{t*) — 2/i G K and 
@ relating to lim^T* w(t*) = J/i € K. 

Equivalently from (|76| ) we deduce that, owing to 
z) > M, M e R, at the least 



lim x 



lim e 3 ^ 4 *)/ 2 . 



Therefore @-© yield for (jj 



lim ^ 

A— »oo // 



0. 



(79) 



(80) 



Along the lines that led to 
satisfies 

lim = 0. 



it is verified that fl75|) 



(81) 



and kO) establish the correctness of (f73]).n 



Eqs. 



6. Ever-expanding quadratic Bianchi 
IX cosmologies 

In this section, we apply Theorems 1, 2 and 3 above 
to some anisotropic quadratic cosmologies and show as 
a corollary of the general theory developed in the pre- 
vious sections that there exist ever-expanding Bianchi 
IX universes in the purely R 2 theory, which may also 
be taken to be regular initially in sharp contrast with 
the situation in general relativity. 

As a warm-up exercise, let us consider first the eas- 
ier cases of models of types I, II and V. 

Case A: V = (Bianchi I, M = 0). 

Theorems 1 and 3 trivially hold but Theorem 2 is 
devoid of any meaning in the present simple case since, 
by setting <q — in all the relevant formulae, without 
loss of generality we obtain 



lim x(t* 

f->T* 



OO, 



T* = - 



InS 

3VCT' 



S 



A 



A 
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Ci = (± ) 2 - fpe 3 * >0; 

/?+ (A) = 2/0-2)0 ln[-*(A)] , 2)0 < 0; 

lim /?+(A) =j/o+ VoT*, 

A — *oo 

/3_(A) = z -iolnH(A)], i < 0, 
lim /3_(A) = zo + zoT*, 

X — >oo 



(82) 



(83) 



where A = d + ^pe 3x ° and S G (0,1)- The case 
when C\ in fl82| ) is not positive leads identically to the 
results of the case C\ > 0. It is probably interesting 
to point out that the case V(y, z) = C = const, which 
is of no physical importance but fulfils the conditions 
for the validity of Theorem 1, yields after a lengthy 
albeit exact calculation by means of elliptic integrals 
the exact finite value of T* with lim t *^T» x(t*) = 00 
and, due to lim t »^ T * v(t ) = Vx G R and 

lime^T* z(t*) = z\ G R. 

Case B: V = exp[4(/3+(A) + \/3/3-(A))] (Bianchi II). 

Eq. @ holds with M = 0, while @ and @) are 
satisfied due to the analytic form of V and to the fact 
that /3+(A) and /3-(A) are monotonically decreasing 
functions of A. Consequently Theorems 1 to 3 are valid. 
Note here that the possibility lim^oo /3+(A) = —00 or 
lim^co /?-(A) = —00 or both is not excluded. 

Case C: V = c 4fj +^ (Bianchi V). 

It is evident that here Theorems 1 to 3 also hold 
and the possibilty lim^oo /3+(A) = —00 may appear. 
Case D: (Bianchi IX). 

We now move on to the Bianchi IX case which from 
the point of view of dynamics is the most interesting. 
The usual Bianchi IX potential is 



V = I e 4/3+(A) 
2 

+ Ip-8/MA) 

+ 4 



- 1 



cosh (4V3/?_(A) 
_ e -2/3+(A) CO sh(2\/3/3_(A) 



Now, on introducing the variables y(t*) and z(t*), by 
virtue of (^) we obtain for the Bianchi IX potential 



W(y,z) 



+ 1 = V'faz), 



(84) 



where V'(y,z) is the positive-definite potential used 
in |l0[ . Therefore M = —3/4, and Theorem 1 retains 
its validity for the Bianchi IX potential. To decide 
whether Theorems 2 and 3 are valid, we need 



V y = 2e 4y 



cosh U*Sz} - 
+ 2e~ 2v c 
V z = 2v / 3[e 42/ sinh(4v / 3z) 



- 2e~*y 

h(2V3z), 

e~ 2v sinh(2\/3z) 



(85) 



In the ensuing developments we consider all the 
possibilities. Unless otherwise stated, all limits are 
taken in the direction y — > — 00, z — 00. 



a) We assume that there exist initial conditions 
x ,x > 0, 2/0,2)0 < 0, z ,z <0 such that 



lim x(t 



and 



= xi G IK 

J^yit*) = -°°> 

lim z(t*) = z\ G R. 

( ._>T* 

Then from (|8q) we obtain 

lim V v = lim(-2c 8M ) = -00, 

Z— >Z\ 

lim V z = lim(-e 2|y| ) = -00. 



(86) 



From (jS6la) it is obvious that (|58| ) is violated for the 
Bianchi IX potential. Furthermore, (|86|a) implies, ow- 
ing to the standard criterion for the divergence of im- 
proper integrals, viz. 



lim (M A e 8H ) = 



A < 1, 



w — > — oc 



that in ( p2|) the integral tends to —00 thereby making 
the square root complex. Hence this case cannot occur. 

b) Suppose now that for appropriate xq , xq > , 
2/o,2)o < 0, z ,z Q < 0, we have lim t »^ T * x(t*) = x± G 
R, lim f ._v T * «(**) = -00, lim t »^ T . y(t*) = yi G R. 
Then (||) yields 



lim V» 



lim(e 



4%/3|zh 



OO, 



2/— >2/i 

lim V z = lim(-e 4v/5|z| ) = -c 



(87) 
(88) 



As in possibility a) above, due to the condition (p8|), 
([59|) ceases to hold for the Bianchi IX potential. More- 
over, since from (|d|) and 



-if 1 / 2 , 



lim i = -if 1/2 (-00), 

i*-»T* 



(89) 
(90) 



we are again led due to (|S8|), so that the square root in 
( po|) becomes complex. Thus this case is also excluded. 

c) Finally, let for some specific initial conditions, 
x ,±o > °. 2/o,2)o < 0, z Q ,z < 0, the result 

lim (x,y,z) — (x%, — 00, — 00), 



be valid. From (85) we deduce 

\imV y = lim(e- 4 ^l+ 4 ^ z l-|-e 2 ^l +2v ^l z l-2e 8 ^l), 



limy, = limV3(-c- 4 ^+ 4 ^l z l+e 2 l«l+ 2 ^ z l) 



(91) 



(92) 



In the following we investigate how ( ]91]) and (|92j) im- 
pinge on the validity of (58) and p9| ) and, conse- 
quently, on the validity of Theorems 2 and 3. 
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There are various possibilities to be considered. 
Firstly, if for certain initial conditions Xo , io > , j/o j 
yo < and zq , io < , we have 

lim(M/l!/l) < V3, 
then (H) and © yield 



limV^ = lim(— exp[8|j/|]) = — c 



(93) 



and 



lim V z = limexp[2|?/| + 2\/3|z|] = oo, 



(94) 

respectively. Eqs. j9^ ) and ( pi|) demonstrate that (5S) 
and (59) are not fulfilled. In addition, by virtue of (92) 
in (|62|), the integral tends to — oo, which renders the 
square root complex. Therefore this case cannot occur. 
Secondly, suppose that for a set of initial conditions, 
Xq , io > , yo, 2/0 < and z , i < , the relation 

hm(|z|/M) > V3 

holds. From (^l|) and (^) we obtain 



lim V y = lim exp 



8|y| + i2|y| 



(95) 



and 

lim V z = lim 



exp 



(2|y|+2V3|i 



= — oo, 



(96) 



respectively. Proceeding as in i) above, we conclude 
that this possibility too is excluded. Lastly, if there 
exist initial conditions, xq , io > , yo , yo < and 
zq, zq < , generating 

]im(\z\/\y\) = v/3, 

© and (H) show that 

lim K = 0. 



lim Vj, = 0, 



(97) 



From (|97|), ( |28| ) and j29| ) we infer by integration in 
the vicinity of T* , as in the proof of Lemma 2, that 
lim t *_>T* y{t*) and lim t »^x* z(t*) are finite, thus ar- 
riving at a contradiction with the supposition made at 
the outset of the present case c). 

Cases a) - c) above immediately show that for 
the Bianchi IX potential and for all initial conditions, 
xq , io > , j/o > Vo < and zq , io < , we have 
limc^y* x(t*) = oo, \imf~>T* y(t*) = j/i 6 1 and 
lim t *^T* = z i £ K. We thus arrive at ( |60| ) 

without having to resort to Lemma 2, and precisely 
such a possibility was announced prior to Lemma 2. 
It is straightforward by following the proof of The- 
orem 3, since now V y and V z are bounded for y S 
[2/1 > 2/o] > z G z o] j to conclude that it holds as ( f73| ) is 
based essentially on lim^^^ x(iu) = lim lu ^ zi x(w) = 
\im t *^T* x(t*) — 00. We thus have 



Corollary 1. For the potentials considered in Theo- 
rems 1 to 3 we obtain for the volume Vl{\) — e 2 ^ x ^ 
by using ( |&i| ) that Vl(A) is a monotonically increasing 
function of A, < A < 00 and, owing to J6?[ ) and (36), 
fulfils 



V L (X = 0) = cxp[2/i(A = 0)] = exp(x + Q, 
lim V L (X) = lim c 2 ' l(A) = 00. 

A — >oo A — >oo 



(98) 



Remark. Had we chosen in (39) the minus sign, we 
would not have been able to prove the basic Lemma 
1. A respective calculation has been carried out and it 
does not lead to any sensible results. Further, the plus 
sign in either (^) or (|d|) or in both would generate 
increasing functions y(t*) or (and) z(t*) , thus prevent- 
ing us from formulating ( pq ) or (|59|) or both in view 
of the analytic form of the Bianchi potentials. Finally, 
it is plausible to use x(t*) as the independent variable 
in (|44|) since it is actually a time variable according to 



(M 



7. Lie point symmetries and Painleve 
analysis of some Bianchi models 

To round off our discussion it is intructive to con- 
sider the question of algebraic integrability of the mod- 
els considered previously. This is certainly far from 
being of purely academic interest since one needs to 
know if the known solutions occupy a large part of 
the phase space or whether there exist regions where 
non-integrability manifests itself in a non-trivial way. 
In such a case usually one has the phenomenon of the 
formation of singular envelopes (cf. g) whereat any 
non-integrable regions are enveloped by the known so- 
lutions. 

One method to determine the integrability of a sys- 
tem is by a performance of the so-called singularity, 
or Painleve, analysis in an effort to examine whether 
or not there exists a Laurent expansion of the solution 
about a movable pole which contains the number of ar- 
bitrary constants necessary for a general solution. Any 
other singularities are not permitted except in the case 
of branch point singularities which give rise to what is 
called the weak Painleve property. A system which is 
integrable in the sense of Painleve has its general so- 
lution analytic except at the pole-like singularity (for 
more details and a nice introduction to the singularity 
analysis of dynamical systems we refer the reader to 
@, ch. 8). 

It is interesting that even for these quadratic mod- 
els where, as we saw earlier, the global dynamics can 
be very different from the corresponding one in gen- 
eral relativity, integrability is not a common feature in 
these examples except for trivial cases. 
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We first examine the system (p7|)-(29) for Lie point 
symmetries. For the computation we use Program LIE 
due to Head g [jj). 

Case A: V(y,z) = (Bianchi I), M = 0. 

The equations are particularly simple, being 



x + a±e 



~3x 







y = 
z = 0. 



(99) 



There are nine Lie point symmetries: dt , d v , d z , td v 
, td z , yd y , zd y , zd y , zd z . We observe that most of 
the symmetry is due to the second and third members 
of (99), both of which are trivially integrable. The first 
has the solution 



x(t) = -~ log •{ {ai/I) 1/2 sinh 



| (2/) 



1/2 



(t-fo) 



where / and to are constants of integration. 
Case B: V = exp[4(y + y/3z)] , (Bianchi II) 
The system is now 



ma 



3x — 24Ke 2x+i( - y+ ^ z ' ) = 0, 
y + i2Kc 2x+i{ y + ^ = 0, 

z + nV3Kc 2x+4 (y+^ = o. 

Upon introduction of the new variables 

v — y + \/3z, w = — v% + z 

the system ( |l00| ) takes the simpler form 

x - 9Ke 3x - 24Ke 2x+4v = 0, 
v + 48K e 2x+4v = 0, 
w = 0. 



(100) 
(101) 



(102) 



The parameter K = 8p, where p was introduced in 
(|l7|). The Lie point symmetries of (102) are dt , d w , 
td w , wd™ . 

We observe a considerable reduction in the number 
of Lie point symmetries compared with Case A. The 



solution of the third of (102) is trivial, but of the first 
two not at all obvious. There is a first integral 



3Ke 6 



YlKc 



2x+4v 



(103) 



Eqs. (102 a) and (102 b) are not in a suitable form for 
applying the Painlcvc test. Under the transformation 



r{T) = 27Ke i ( 2x+v '> 
s(T) = 27Ke 3x , 



T 



Ut 



(104) 



(the constants are chosen for later numerical conve- 
nience and do not affect the essence of the analysis) 
we obtain the system 



tt — r 2 + 6r 2 s = 0, 



s d + 6r = 0. 



(105) 



(We continue to use an overdot to indicate differentia- 
tion with respect to the new time variable, T .) Making 
the usual substitution 



OLT 1 



(106) 



in (105), we obtain p = —6,q = —2 and p > —6,q = 
—2 as the possible exponents. For the former the lead- 
ing order behaviour is 



.I r -6 

2 ' ' 



(107) 



When we examine the system (105) for resonances as- 
sociated with (107) we find that they are -1, 2 and 
5 (1 ± i-v/35) ■ The system ( |105| ) will not be integrable 
in the sense of Painleve and, although the transforma- 



tion (104) is not homeographic, we can infer the same 
for the original pair of equations. The meagre num- 
ber of Lie point symmetries associated with the pair of 
equations, just G\ , reinforces this conclusion. 

Case C: V = e iy (Bianchi V). 

The system is now 



9Kc 3x - 24Kc 2x+4y = 0, 
y + 12Kc 2x+iy = 0, 
z = 0. 



(108) 



The symmetries are just those in Case B above with w 
replaced with z . The integration of the third equation 
is trivial and the first two have the first integral 



J 



2 x 



2y 2 - 3Ke 3x - l2Ke 2x+4y . 



(109) 



As was the case with (102), (|109| ) is not in a form 
suitable for the Painleve analysis. We introduce the 
transformation 



r(T) 
s(T) 



_ _ e 2x+4y 



T = i(l/9K) 1/2 t, (110) 
in which the numerical coefficients are chosen for arith- 



metrical convenience. The system (109) becomes 



r 2 + 3r 3 + 5r 2 s = 0, 
ss - s 2 + 2r 2 = 0. 



(Ill) 



Both indices are —2 and the only possible leading order 
behaviour requires that all terms be dominant. The 
resonances are at — 1 , 2 and | (1 ± \/l9) , and we 
conclude that this system is also not integrable. 
Case D: Bianchi IX. 

The only symmetry which this system possesses is 
dt , that is, the system is invariant under time trans- 
lation. It is not integrable. It is known that the non- 
integrability of this model in general relativity associ- 
ated with the BKL oscillatory behaviour of the scale 
factors in the contracting direction disappears in the 
framework of the quadratic theory derived from the La- 
grangian R + aR 2 and in all f(R) Lagrangian gravity 
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theories that contain an Einstein term. This is due to 
the conformal equivalence of such theories to general 
relativity with a scalar field matter source and this case 
is known to have a monotonic approach to the initial 
Bianchi IX singularity ||. Here, however, we see the 
interesting result that in the purely quadratic theory 
(that is, without the Einstein term) some sort of non- 
integrability returns. It remains to be seen if this lack 
of integrability of the Bianchi IX spacetime found here 
is related to a possible chaotic behaviour in the space- 
time geometry similar to that which occurs in general 
relativity. If true, such a result would imply, for in- 
stance, that the pure R 2 theory is in a sense closer to 
general relativity than other quadratic extensions of it 
containing the linear Einstein term. 

8. Conclusion 

The behaviour of vacuum quadratic cosmologies in 
the expanding direction through the Hamiltonian ap- 
proach advanced here can be generalized to include 
matter fields, and this could be the next problem 
to tackle in this respect. How does the inclusion of 
a scalar field or a perfect fluid affect the isotropiza- 
tion theorem proved in Sec. 5? Is the non-recollapse 
of the Bianchi IX model in quadratic gravity stable 
to matter inclusion or to perturbations of the purely 
quadratic scalar curvature Lagrangian? In general 
gelativity all homogeneous, vacuum or matter-filled, 
closed cosmologies satisfying the usual energy condi- 
tions recollapse (cf. ]l0||, pdj]). We have shown that 
the purely quadratic diagonal Bianchi IX models are 
ever-expanding (i.e., volume-increasing) and so they 
provide a counterexample to a possible closed universe 
recollapse conjecture in higher-order gravity theories. 

Another question which can be of interest is to anal- 
yse these systems in terms of the expansion- normalized 
variables of Wainwright (cf. 0) and compare the re- 
sults with those arrived at here. This would nicely 
complement the isotropization theorem proved in this 
paper about the scalar shear variable. 
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